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ABSTRACT 



We present some important conclusions from recent calculations pertaining to the collapse of rotating molecular cloud cores with axial 
symmetry, corresponding to evolution of young stellar objects through classes and begin of class I. Three main issues have been 
addressed: (1) The typical timescale for building up a preplanetary disk - once more it turned out that it is of the order of one free-fall 
time which is decisively shorter than the widely assumed timescale related to the so-called "inside-out collapse"; (2) Redistribution of 
angular momentum and the accompanying dissipation of kinetic (rotational) energy - together these processes govern the mechanical 
and thermal evolution of the protostellar core to a large extent; (3) The origin of calcium-aluminium-rich inclusions (CAIs) - due to 
the specific pattern of the accretion flow, material that has undergone substantial chemical and mineralogical modifications in the hot 
(k. 900 K) interior of the protostellar core may have a good chance to be advectively transported outward into the cooler remote parts 
4 AU, say) of the growing disk and to survive there until it is incorporated into a meteoritic body. 
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1. Introduction 

Stars are thought to form from regions within molecular clouds 
that for some reason became gravitationally unstable and started 
to collapse under the influence of their own gravitation al at- 
traction (see, e.g.. lLarsonll2003tlMcKee & Ostriketll2007l). It is 
clear that angular momentum (e . g. jBod enheime rj 1 1 9951; iTo hline . 
2002) and magnetic fields (e.g., Pudri tz et all 120081) play an im- 
portant role during many stages of protostellar collapse. Angular 
momentum is particular important because it is responsible for 
the formation of accretion disks that are the birthplaces of plan- 
etary systems. Therefore, the study of rotating collapse is in- 
evitable if one intends to understand the formation of the Solar 
System and other planetary systems. 

The earliest such studies are a number of analytic a pproaches 
for rotating collap se without dTerebev et all 1 19841) and with 
dAllen et all 120031) magnetic fields and analytic studies of disk 
formation (see Cassen, 1994, and references therein) and sim- 
ple models for build-up and evolution of accretion disks (e.g., 
iLin & Pring5ll990llNakamoto & Nakagawalll994l) . These suf- 
fer from that they introduce a number of approximations on the 
nature of the collapse process that appear plausible but lack jus- 
tification. 

In the past two decades various attempts have been made, 
therefore, to follow the 2-D and 3-D collapse of rotating pro- 
tostellar fragmen ts by means of extended nume r ical calcula- 
tions (e.g.. iBossl [19891; iBoden heimer et all 1 19901; lYorke et all 
1993t lSaigoet7lT T2008. to give a few examples) based on grid- 
methods, by a solution meth od based on a seri es development 
into orthogonal polynomials dTscharnuteilll987|). and by studies 
based on the SPH-method (e.g., Walc h et al. , 2009). However, 
because of lack of spatial resolution almost all of these studies 
followed only the first collapse phase up to the formation of the 



first core. Neither the second collapse of the core to stellar di- 
mensions (^ 1 Rq) could be followed, nor the formation of the 
inner par t of the accretion di sk on size scales <k 1 A U. Only the 
model of iTscharnutetj d 19871) and the recent study of Saig o et ail 
(2008[) ac hieved sufficient re solution to follow the second col- 
lapse, but lSaigo et al.l (120081) neglects most of the basic physics 
and considers only an over-simplified equation of state. 

Many of the investigations mainly focus on the individual 
stability behaviour of the rapidly spinning, flattened, quasi-equi- 
librium "cores" that form in the innermost optically thick parts 
of collapse flows. It has long been known that, for a rotating self- 
gravitating equilibrium configuration, the ratio of rotational over 
gravitational energy, 9 :- E TOt / LEgrav Jjis the appropriate param- 
eter to indicate the onset of secular, for 6 0.14, and dynamical, 
for 6 =i 0.27, instability. This is a standard result for homoge- 
neous equilibria, but interestingly enough, in a series of papers 
Bode nheimer & OstrikeA Il970ill973h lOstriker & Bodenheimeti 
1973 were able to show that these two critical numbers still hold 
for more general rotating polytropes, almost independently of 
the poly tropic index, n. 

Due to this remarkable finding extensive 3-D calculations 
have been conducted in order to explore the fate of rapidly 
spinning "naked" protostellar cores by following the onset of 
the dynamical instability and its development into the nonlin- 
ear regime. Much effort has gone into clarifying the growth 
and saturation of the self-gravitating non-axisymmetric modes 
(bars, spirals, mixed forms) and their mutual interactions which 
drive a ver y effective redistribution of the cores' angular momen- 
tum ( e. g.JPickett et aUll996h iToman et all ll998tlPickett et all 
ll997Ulmamura et alll2000l 

In this letter we present some results of a new model calcu- 
lation that follows for the first time simultaneously the evolution 
of a rotating Bonnor-Ebert-sphere-like initial state with ~ 1 Mq 
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1 Usually written fj : = T/ \ W | , but we wish to avoid confusion with T , 
the temperature, and /3, the scaling parameter of the turbulent viscosity. 
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Table 1. Starting parameters. 



Quantity Symbol Value Dimension 



mass 




M 




1.037 


M© 


angular momentum 




L 




2.682- 10 53 


2 -1 

g cm - s 


radius 




R 




1.200-10 17 


cm 


spin angul. veloc. 




Q. 




3.160-10~ 14 


s- 1 


centrifugal barrier 






1.504-10 15 


cm 


mean density 




P 




2.850- 10- 19 


gem -3 


mean free-fall time 


is 




3ir 

nap 


1.247- 10 s 


yr 


central density 




Pc 




5.326-10- 18 


gem" 3 


central free-fall time 


fff.c 


-V 


371 

32Gp c 


2.885- 10 4 


yr 


temperature 




T 




10 


K 


ratio rot./grav. energy 


6 


__ £roi 
~ lEgravl 


2.437- 10- 3 





initial mass through axially-symmetric first and second collapse 
to nearly stellar central densities and the early build-up and evo- 
lution of the associated accretion disk. The maximum resolu- 
tion achieved amounts to about 0.03 Rq in the central parts. The 
adaptive grid is chosen in a way that the innermost rigidly rotat- 
ing homogeneous sphere always contains a fraction of 10~ 6 of 
the total mass. The simulation covers the evolutionary phases of 
young stellar objects corresponding to class and early class I. 
It clearly demonstrates the rapid co-formation of a compact stel- 
lar object and a very extended accretion disk within a period 
of no more than ~ 1.1 free-fall times. The model includes all 
of the essential physics of the problem: hydrodynamics, gravity, 
radiative transfer (in Eddington approximation), and a realistic 
equation of state including effects of sublimation of ice and dust 
on the opacity. 

By using a fully implicit numerical method (see Section [2]l 
we are now able to follow the collapse of a slowly rotating, 
slightly Jeans-unstable molecular cloud fragment with axial 
symmetry. 

2. Methods 

The most frequently applied nu merical meth ods - such as nested 
grids by Saigo et al. or SPH bv lWalch et al.l - have in common 
that they are of explicit type ("forward" time differences), i. e., 
they are subjected to the Courant-Friedrichs-Lewy (CFL) con- 
dition for the timestep so as to warrant numerical stability. In 
fact, the CFL-timestep becomes prohibitively low with respect 
to the accretion time, if the structure of the hot stellar core is to 
be sufficiently well resolved. This is why explicit methods of any 
kind are genuinely unsuited for attacking the overall problem of 
stellar formation in a consistent way. 

Unfortunately, appropriate implicit methods ("backward" 
time differences) have been worked out only for 1-D and 2-D 
collapse problem s with spherical and axial symmetry, respec- 
tively (see, e.g., IWuchterl & Tscharnuterl 120031; iTscharnuterl 
1987). Although implicit numerical schemes do not suffer from 
the CFL-condition, they demand, after all, the solution of a huge 
system of non-linear equations which is computationally expen- 
sive. Thus, implicit methods are useful only if the physical pro- 
cesses to be considered exhibit a hierarchy of timescales and the 
system as a whole evolves into a quasi-stationary state. If the 
longest timescale exceeds the smallest one by a large amount, 
the use of implicit methods is mandatory, as is in our case, where 
the accretion time is several orders of magnitude longer than the 
oscillation period of the stellar core. 
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Fig. 1. Meridional cross sections. Displayed are equi-density 
contours, the spacing being 0.5 dex, and the velocity field. Red 
"spots" indicate expansion (cf. Fig. [3^). The closed heavy line 
marks the accretion shock, the arrow at the lower right corner 
of the four panels represents the respective maximum infall ve- 
locities. Upper row: snapshots slightly after B and before C, re- 
spectively. Lower row: snapshots for D and E (cf. Fig. |2). At E 
(right panel) two shocks, an outer and an inner one, bounding the 
growing "pre-planetary" disk and the "stellar" core, respectively, 
coexist (not displayed here, but cf. Fig.©. 



For a detailed description of the method which t he implicit 
2-D h ydrodynamical code is based on we refer to Tsch arnuterl 
(1987). For the equation of state, dissociation of H2 and ioni- 
sation of H and He is taken into account; the opacities at low 
temperatures are dominated by dust with and without ice man- 
tles. In a spherical polar coordinate system, we discretize the 
variables on a staggered radial grid (256 gridpoints) and repre- 
sent the dependence on the polar angle by choosing a Legendre 
expansion (up to 27 coefficients) for each primary variable. The 
discretized equations are written i n conservation f orm on a self- 
adaptive radial grid according to Dorfi & Druryj d 1987b . shock 
fronts are smoothed out by artificial (tensor-) viscosity. 

Concerning the coefficient of turbulent viscosity, v tul , 
we have adopted the so-called y8-viscosity prescription of 
iDuschl et al. (2000), an extension of the well-known, but like- 
wise heuristic, a-viscosity for application to self-gravitating 
disks. To be specific, v tul := j6r 2 Q, with r being the radial co- 
ordinate, Q. the (local) angular velocity, and f3 = 10~ 4 . . . 10~ 2 
the inverse of the critical Reynolds number indicating the onset 
of turbulence. In our calculation we have chosen (3 - 10~ 2 . 



3. Results and conclusions 

3. 1 . Initial configuration and the first collapse phase 

Table [T] lists a set of appropriate starting parameters leading to 
the collapse of a rotating protostellar cloud fragment. The (op- 
tically thin) Bonnor-Ebert-sphere-like initial configuration ex- 
hibits already a moderate density concentration toward the cen- 
tre and is assumed to rotate like a rigid body. The revolution 
period, 2n/Q. = 6.3 ■ 10 6 yr, adopted in our calculation is large 
compared to the mean free-fall time. Hence, the centrifugal bar- 
rier is situated at a distance from the centre of only about 1 % of 
the cloud radius, and the collapse flow develops almost perfectly 
with spherical symmetry for most parts of the cloud. 
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Fig. 2. Evolution of the core, (a): central density vs. central tem- 
perature, (b) ratio of rotational over gravitational energy vs. core 
mass, (c) core mass vs. time. For the meaning of labels A through 
E see SectionP 



After about one mean free-fall time the temperature rises in 
the centre of the fragment and a quasi-hydrostatic core is going 
to form. It is thus convenient to reset the clock to zero when 
the optical depth of the collapsing fragment, counted from the 
outer edge to the centre, exce eds 2/3 for the first time (e. g., 
IWuchterl & Tscharnuterll2003l) . This event may be attributed to 
the beginning of the protostellar evolution proper. All ages sub- 
sequently given are relative to this instant. 

Shortly after age zero an accretion shock forms, that marks 
the natural boundary of a "core", a flattened quasi-hydrostatic 
rotating structure that is pressure-supported parallel and, in 
essence, centrifugally-supported perpendicular to the axis of ro- 
tation. 

As the "core" we shall refer to the innermost subsonic re- 
gion of the collapsing fragment, regardless of whether or not the 
accretion shock is present. During the main accretion phase the 
flow becomes indeed highly supersonic and a strong shock is 
present at any polar angle (cf. Fig. [T]). However, the retarding 
effect of the centrifugal forces significantly reduces the veloc- 
ity components perpendicular to the axis of rotation. Particularly 
close to the equatorial plane, this effect will sooner or later in the 
evolution weaken the local shock strength considerably or even 
make the flow discontinuity disappear completely. This will hap- 
pen within only a few initial mean free-fall times. Then infall of 
matter from the envelope is eventually going to cease, whereas 




Fig. 3. Evolution of the core from its quasi-stationary stage to 
the end of the second collapse. Spatial distribution of (a) the ra- 
dial velocity component, u r , (b) the temperature, T, in the equa- 
torial plane for three instants of time indicated by dash-dotted 
(2301 yr), dashed (2948 yr), and solid lines (3070yr), respec- 
tively; crosses mark the gridpoints. 



the "core" gradually turns into a veritable protoplanetary accre- 
tion disk. 

The current simulations do not only cover the collapse 
proper, but also the formation and the beginning of the ensuing 
accretion phase of the disk-like core. 

3.2. Formation and growth of the disk-like core 

Figure [2] shows (a) the density-temperature-diagram for the very 
centre of the collapsing cloud, (b) the stability parameter 6 C as a 
function of the relative core mass M c /M lot , (c) the time depen- 
dence of the relative core mass. The labels A through E mark 
characteristic stages of the evolution: 

A: The core starts to grow in mass; 

B: continuation of the accretion after thermal relaxation caused 
by opacity effects (sublimation of ice mantels); 

C: sublimation of the dust grains commences off-centre near the 
axis of rotation as a "hot polar cap", leaving a growing opac- 
ity gap and violent thermal relaxation effects; 

D: start of the second collapse due to dissociation of H2; 

E: formation of the "stellar" core. 

As the first important result, we have found a typical rise 
time, M c /M c , of the core's mass to be only a small fraction 
(1-2%) of the mean free-fall time, i.e., only a few thousand 
years. Figure [2J: shows an almost constant accretion rate of 
3 • 1CT 5 Mq yr _I lasting for about 2 600 yr, from slightly after "A" 
until "D" where the second collapse sets in. The dots between 
"D" and "E" represent the individual models tracing the dynam- 
ical transition (the "second collapse") from the disk-forming first 
core to the second "stellar" core. 

Figure Q] demonstrates the increasing geometrical dimen- 
sions and the changing shape of the first core, which occurs 
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Fig. 4. Meridional cross sections (cf. Fig. [TJ. Formation of the 
"stellar" core on a very short timescale. Plotted are three snap- 
shots of the final shock evolution at (from above) to - 10 d, to, 
and to + 7 d, respectively, with fo = 3070 yr. Note that now the 
geometrical dimensions are only several solar radii (Rq). 



within a few thousand years; the bulge in the polar direction, in- 
dicated by the route of the accretion shock in the vicinity of the 
rotational axis, is the result of a rather intricate interplay between 
redistribution of angular momentum, the accompanying genera- 
tion of entropy (viz. heat), and transport of energy. The net effect 
are rising temperatures and lower rotation rates in the inner parts 
of the core, which is necessary to trigger the second collapse. As 
a matter of fact, test calculations have shown that without a cer- 
tain amount of angular momentum transport the core would stay 
too cool and, hence, become prone to the onset of destructive 
gravitational instabilities that trigger binary formation. 

A glimpse to Fig. [2p. reveals that the first core cannot escape 
from running into the regime of secular instability. The ratio, 
6 - E mt /\E grdv \, quickly rises above the critical number of 0.14, 
so that already a small amount of dissipation will lead to sym- 
metry breaking: The originally axisymmetric core will take on a 
triaxial shape and presumably evolve into a distinctive bar/spiral 
configuration exerting gravitational torques which, in turn, re- 
sults in an enhanced redistribution of angular momentum. This 
finding may serve as motivation for our rather efficient turbulent 
(jS-)viscosity with /3 = 0.01. 



3.3. Formation of the "stellar" core 

The second collapse is a result of the interplay between ther- 
modynamics, redistribution of angular momentum as a dissipa- 
tive process, and energy transport. Their combined effect is il- 
lustrated by Fig.[T] particularly in the central parts, angular mo- 
mentum transport tends to defuse extrem flattening, while the 
heat input by the accompanying dissipation of rotational energy 
creates increasing pressure forces. As a consequence, the pres- 
sure distribution is becoming more spherical, and after some "in- 
cubation" period of slow contraction (between "C" and "D" in 
Fig. |2]) dissociation of H2 eventually leads to dynamical collapse 
("D"-"E"). Figure I3J) shows that the stability parameter, 8 C , ap- 
proaches the critical number of 0.27 indicating dynamical insta- 
bility. Interestingly enough, our axisymmetric model suggests 
the formation of a self-gravitating doughnut-like structure of the 
density distribution (cf. Fig.|4|. 

3.4. An intermediate hot solar nebula? 

The model provides for the first time the initial temperature and 
density structure of the accretion disk and the flow pattern in 
the disk and its surroundings, including the accretion shocks, on 
scales down to much less than 1 AU. Figure [3}? shows a rather 
hot inner disk region extending out to about 4 AU. This structure 
is not seen in earlier calculations because of lack of resolution; it 
lasts for a couple of 1 000 yr. This inner portion is hot enough for 
formation of materials like those found as calcium-aluminium 
rich inclusions (CAIs) in meteorites; it might be the "hot solar 
nebula" that cosmochemists always advocated for. 

From the onset of disk formation on there is both accre- 
tion and a large-scale transport of matter (with velocities 10- 
50ms~') close to the disk's midplane from the hot part out- 
ward to distances of several AU from the centre (cf. Figs. Q] 
and [3^). Material from the hot region may be mixed across the 
inner ~ 5 AU and part of this material may survive till the onset 
of planetesimal formation, since outflow close to the midplane 
conti nues to exist in later p hases, though with reduced velocity 
(e.g., iKeller & Gail l2004f) . The model therefore seems to of- 
fer an explanat i on ho w the short, maybe as short as 20000yr 
dJacobsen et al.L 120081) . CAI-forming period in the Solar Nebula 
is related to the earliest evolutionary phase of the accretion disk. 
If true, CAI formation follows immediately the second collapse 
and is a valid indicator for the formation time of the protoplane- 
tary disk. 

Acknowledgements. This work has been supported by the Forschergruppe 759, 
"The Formation of Planets: The Critical First Growth Phase" of the Deutsche 
Forschungsgemeinschaft (DFG). 



References 

Allen, A., Li, Z.-Y., & Shu, F. H. 2003, ApJ, 599, 363 

Bodenheimer, P. 1995, ARA&A, 33, 199 

Bodenheimer, P. & Ostriker, J. P. 1970, ApJ, 161, 1101 

Bodenheimer, P. & Ostriker, J. P. 1973, ApJ, 180, 159 

Bodenheimer, P., Yorke, H. W., Rozyczka, M., & etal. 1990, ApJ, 355, 651 

Boss, A. P. 1989, ApJ, 345, 554 

Cassen, P. 1994, Icarus, 112, 405 

Dorfi, E. A. & Drury, L. O. 1987, J. Comp. Phys., 69, 175 

Duschl, W. J., Strittmatter, P. A., & Biermann, P. L. 2000, A&A, 357, 1123 

Imamura, J. N., Durisen, R. H., & Pickett, B. K. 2000, ApJ, 528, 946 

Jacobsen, B., Yin, Q., Moynier, F, et al. 2008, Earth Plan. Sci. Lett. 272, 353 

Keller, C. & Gail, H.-P. 2004, A&A, 415, 1177 

Larson, R. B. 2003, Reports on Progress in Physics, 66, 1651 

Lin, D. N. C. & Pringle, J. E. 1990, ApJ, 358, 515 

McKee, C. F. & Ostriker, E. C. 2007, ARA&A, 45, 565 

Nakamoto, T. & Nakagawa, Y. 1994, ApJ, 421, 640 



Werner M. Tscharnuter et al.: Protostellar collapse: rotation and disk formation 



5 



Ostriker, J. P. & Bodenheimer, P. 1973, ApJ, 180, 171 
Pickett, B. K., Durisen, R. H., & Davis, G. A. 1996, ApJ, 458, 714 
Pickett, B. K., Durisen, R. H., & Link, R. 1997, Icarus, 126, 243 
Pudritz, R. E., Banerjee, R., & Ouyed, R. 2008, ArXiv e-prints 
Saigo, K., Tomisaka, K., & Matsumoto, T. 2008, ApJ, 674, 997 
Terebey, S., Shu, F. H., & Cassen, P. 1984, ApJ, 286, 529 
Tohline, J. E. 2002, ARA&A, 40, 349 

Toman, J., Imamura, J. N., & Pickett, B. K. 1998, ApJ, 497, 370 
Tscharnuter, W. M. 1987, A&A, 188 

Walch, S., Burkert, A., Naab, T, & Gritschneder, M. 2009, ArXiv e-prints 

Wuchterl, G. & Tscharnuter, W. M. 2003, A&A, 398, 1081 

Yorke, H. W., Bodenheimer, P., & Laughlin, G. 1993, ApJ, 411, 274 



